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ABSTRACT 


Volatility in stock markets is a routine phenomenon. It has been an interesting and challenging area as well for the Stock 
Market investors since long. A lot of researches were carried out in the past taking help of various mathematical models 
but still it is an ongoing process. This study was undertaken to look into the various mathematical models applied to 
analysis the volatility in stock markets. Various equations, tests for stationarity, statistical tools and GARCH family models 


for heteroscedasticity were undertaken in the study. 
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INTRODUCTION 


Some people believe that stock market is a risky market to invest their hard earned money. On the other hand, others want 
to see the stock market bear and bull and to be part of that hassle. Now, this is all volatility in the market which makes this 
place unpredictable and risky one. It is believed that volatility and return are intermingled and volatility has an impact on 
returns of the stock market. Mittal (2009) analyzed volatility & return of the capital market of India from the perspective of 
understanding market behavior. It was observed that annual volatility declined from 2000 to 2007. Since then, however 
there was sharp increase in volatility. It is seen how the market responds to different happenings either positive or negative 
and how it effects the returns of the market. The BSE & NSE are two major stock exchanges in India as most of the share 
transactions are done by the investors in these two exchanges. Derivatives were introduced in the Indian stock market in 
2001to control the volatility of Indian stock market. Bhatia & Jindal (2020) described various types & dimensions of 
volatility in stock markets. Nisha (2014) applied various mathematical & statistical models for analyzing stock market 


return and volatility in her study. 
Statistical Tools used to Analysis the Trends in Stock Market Returns and Volatility Patterns 


To calculate the returns, logarithmic difference between two periods is primarily taken by applying the following: 


R,Uin P, — In P,_,) * 100 ; 


Where Rt is the return in period t, Pt and Pt-1 are the daily closing prices of the index at time t and t-1 


respectively. 
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Unit Root Test 


For testing stationarity, let us consider an AR (1) model: 


Vi = Pi¥iat & i 


The simple AR (1) model indicated in above equation is classified a random walk model. In this AR (1) model if 
Ipil<1, then the series is I(0) i.e. stationary in level, but if p= 1 then there exist what is called unit root problem. In other 
words, series is non-stationary. Most economists think that differencing is warranted if estimated p > 0.9; some would 


difference when estimated p > 0.8. Besides this, there are few formal means of testing for stationarity of a series. 
Augmented Dickey Fuller Test 


Dickey-Totaler test involve estimating regression equation and carrying out the hypothesis test. The simplest approach to 


testing for a unit root is with an AR(1) model. AR(1) process: 


Fe = C+ P¥i-at ©: 3 


Where c & p are parameters; is assumed to be white noise. If -l1<p<1, then y is a stationary series while if p=1 , y 
is a non-stationary series. If the absolute value of p is greater than one, the series is explosive. Therefore, the hypothesis of 
a stationary series involves whether the absolute value of p is strictly less than one. The test is done by estimating an 


equation with y,.;subtracted from both sides of the equation: 


Ay, =€+ ¥¥4+ & 4 


The DF test is valid only if the series is an AR(1) process; if the series is correlated at higher order lags, 
assumption of white noise disturbances is violated. The ADF controls for higher order correlation by the adding lagged 


difference terms of dependent variable to right-hand side of regression: 
AY, = C+ ¥¥p-1 + Gy AVy-1 + G2A¥, 2+ ~.. + GAY,» + &; 
This augmented specification is then tested for in this regression. 
Ay: yr =0 
A acF¥ = 0 
The empirical research work done earlier has used this tool in their work (Kaur, 2004; Abdalla, 2012; Bordoloi& 
Shankar, 2008 and Karmakar, 2007). 


Phillips—Perron Test 


In statistics, the Phillips—Perron test (named after Peter C. B. Phillips and Pierre Perron) is a unit root test. That is, it is used 
in time series analysis to test the null hypothesis that a time series is integrated of order 1. It builds on the DF test of the 


null hypothesis é = Oin 
Ayt = 5y,_1 + pat 6 
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Where A is the first difference operator. Like the augmented Dickey—Totaler test, the Phillips—Perron test 


addresses the issue which process generating data for ¥y might have a higher order of auto correlation than is admitted in 
test equation — making Wy_4 endogenous & thus invalidating the Dickey—Fuller t-test. Whilst the augmented Dickey— 


Fuller test addresses this issue by introducing lags of AV; as regressors in test equation, the Phillips—Perron test makes a 


non-parametric correction to the t-test statistic. The test is robust with respect to the unspecified autocorrelation & 
heteroscedasticity in the disturbance process of test equation. (Kaur, 2004; Bordoloi& Shankar, 2008 and Karmakar, 2007) 


used the same for their research work. 
Kwiatkowski-Phillips-Schmidt-Shin (KPSS) Test 


In econometrics, Kwiatkowski—Phillips-Schmidt-Shin (KPSS) tests are used for testing a null hypothesis that an 





observable time series is stationary around a deterministic trend. The series is expressed as the sum of the deterministic 
trend, random walk & stationary error, the test is the Lagrange multiplier test of the hypothesis that the random walk has 
zero variance. KPSS type tests are intended to complement unit root tests, such as the Dickey—Fuller tests. By testing both 
the unit root hypothesis and stationarity hypothesis, one can distinguish series that appear to be stationary, series that 
appear to have a unit root & series for which the data (or the tests) are not sufficiently informative to be sure whether they 


are stationary or integrated. (Gupta and Basu, 2007), earlier literature mentioned the use of KPSS. 
Autocorrelations and ACF (k) 


Autocorrelation is one of statistical tools used for measuring the dependence of the successive terms in a given time series. 
Hence, it has been widely used to measure dependence in successive share price changes. Autocorrelation has been the 
basic tool used to test the weak form of EMH. The autocorrelation function ACF(k) for time series Yt and the k lagged 


series Yt-k is defined as: 


ACF(K) = vpea—lFt— sent Tt 


rn ,0,- v2 ; 
Where Y is considered as overall mean of series with n observations. The SE of ACF (k) is given by: 
S@acriK) = = 
in—kh 8 
When n is sufficiently large (n >50), the approximate value of standard error of ACF (k) is given by: 
S@acriK) = = P 
To test whether ACF (k) is significantly different from zero, the distribution of t is used; 
ACF(K) 
5 acr(K) 10 


As it is true for random walks, trends are also characterized by the extremely high autocorrelation. For both 
random walk series & series with trends, autocorrelation ACF (k) are very high & decline slowly as the lag value (k) 


increases. At the same time the ACF (k) of the first difference series (price changes or returns) are statistically insignificant 
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when the series is a random walk series. A random walk series drifts up and down over time. In some situations it may be 
difficult to judge whether a trend or drift is occurring. Hence to determine whether a series has a significant trend or 


whether it is a random walk, the t-test is applied on series of first differences. (Bordoloi & Shankar, 2008). 
Heteroscedasticity 


One of most important issues before applying Generalized Autoregressive Conditional Heteroscedasticity (GARCH) 
methodology is to first examine residuals for evidence of heteroscedasticity. To test for the presence of heteroscedasticity 
in residuals of KSE index return series, the Lagrange Multiplier (LM) test for ARCH effects proposed by Engle (1982) is 
applied. In summary, test procedure is performed by first obtaining the residuals t e from ordinary least squares regression 
of conditional mean equation which might be an autoregressive (AR) process, moving average (MA) process or a 
combination of AR & MA processes; (ARMA) process. For example, in ARMA (1,1) process the conditional mean 


equation will be as: 
1, OF +E + Oey u 
After obtaining the residuals ¢ e , the next step is regressing the squared residuals on a constant and q lags as: 


€f = Ag+ ayer, + aze; 2+ ..+ a,eF 


Le 


12 


The null hypothesis that there is no ARCH effect up to order q can be formulated as: 


Hy: ay _ ~_ %a_9 


Against the alternative: 


Ay. a@;>0 


For at least one i= 1, 2, ...,q 


The test statistic for joint significance of q-lagged squared residuals is the number of observations times the R- 
squared (TR?) from regression. TR’is evaluated against y7(/(q) distribution. This is asymptotically locally most powerful 


test. 
The Generalized Autoregressive Conditional Heteroscedasticity (GARCH) Model 


The conditional variance is represented as a linear function of its own lags in this model. The simplest model specification 


is the GARCH (1,1) model: 


Mean Equation r,= pt &, B 


Variance Equation a7_w+ a, &%_,+ By ory 14 


where a > Qeanda, >Oand § > O.and 
% =return of the asset at timet 


B= average recurn 
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£, = residual returns, defined as: 


ey = ore; 


Where z, is standardized residual returns (i.e. iid random variable with zero mean & variance 1), and oO, is 
conditional variance. For GARCH (1,1), the constraints a sand B;> 9 are needed to ensure oO, is strictly positive. In this 
model, the mean equation is written as a function of constant with an error term. Since a”, is the one —period ahead forecast 
variance based on past information, it is called conditional variance. The conditional variance equation specified as a 


function of 3 terms: 
e Aconstant term: @ 


e News about volatility from the previous period, measured as the lag of the squared residual from mean equation: 


es 
t-1 ( ARCH term) 


e = Last period forecast variance: Te ( GARCH term) 


The conditional variance equation models time varying nature of volatility of residuals generated from mean 
equation. This specification is often interpreted in a financial context, where an agent or trader predicts this period’s 
variance by forming a weighted average of a long term average (constant), the forecast variance from last period (GARCH 
term), and information about volatility observed in previous period (the ARCH term). If the asset return was unexpectedly 
large in either the upward or the downward direction, then the trader will increase the estimate of the variance for the next 
period. The general specification of GARCH is, GARCH (p, q) is as: 


ot =@+ dja Fy 4+ Di i%s 4 e 


Where, p is number of lagged a’ terms and q is the number of lagged €” term 
The Exponential GARCH (E-GARCH) Model 
This model captures asymmetric responses of time-varying variance to shocks & at same time, ensures that variance is 


always positive. It was developed with the following specification: 


era jo ee 


Ln(o7) = w+ B,Ln(o7_,)+ cr | 








Or-1 Tr r-4 


16 


Where y is asymmetric response parameter or leverage parameter. The sign of y is expected to be positive in most 
empirical cases so that a negative shock increases future volatility or uncertainty while a positive shock eases the effect on 
future uncertainty. In macroeconomic analysis, financial markets & corporate finance, a negative shock usually implies bad 
news, leading to a more uncertain future. Consequently, shareholders would require a higher expected return to compensate 
for bearing increased risk in their investment. Above Equation is an E-GARCH (1,1) model. Higher order E-GARCH 
models can be specified in a similar way; E-GARCH (p, q) is as follows: 


i} ss f2@ Fee 
‘si £ Ve bal ed 
17 











Ln(a?) = wi? =i B; Ln( oF =) + ac; } |= 


a. 


www.iaset.us editor @iaset.us 


34 Vinod Bhatia & Ritu Goel Jindal 


The Threshold GARCH (T-GARCH) Model 
Another volatility model commonly used to handle leverage effects is the threshold GARCH (or T-GARCH) model. In T- 
GARCH (1,1) version of model, the specification of conditional variance is: 


oy Z Z a 
Op _ Wt Oye yt Veep + Bi oe 18 


Whats d-_4 is adummy variable, that is: 


d,. = { 1ife:-1 50 adnews } 
-1 Oif =a =p, good news 

The coefficient y is known as the asymmetry or leverage term. When y = 0, the model collapses to the standard 
GARCH forms. Otherwise, when the shock is positive (i.e., good news) the effect on volatility is 0), but when news is 
negative (i.e., bad news) effect on volatility is a),y. Hence, if y is significant and positive, negative shocks have a larger 
effect on 6”, than positive shocks. This way, T-GARCH (p,q), the conditional variance equation is specified as follows: 

= Y ri P ra 
oO; = @t De hh + Fide 1 )e4 + Yi-4 Boy; 19 

aj, ¥i and6; arenont—negativep arameterssatisfyingconditionssimilartathoseaof GARCHmodel., 


(Akgiray, 1989; Ballie and DeGennaro, 1990; Lamoureux and Lastrapes, 1990; Corhay and Tourani, 1994; Geyer, 1994; 
and Sakata and White, 1998) used ARCH, GARCH and their extension models in their research work. 


CONCLUSIONS 


The various mathematical, statistical & econometrics models considered in the study & depict the wide picture of 
analyzing the volatility of stock market. It included simple statistical methods those are used to forecast the volatility. 
Another important segment is of testing stationarity of the time-series data for which Unit Root Test, ADF Test and PP 
Test are applied. Econometric model i.e. KPSS Test had been very popular for testing a null hypothesis in the process of 
checking stationarity of data. Autocorrelation had been another important tool for measuring the dependence of successive 
terms in a given time-series proceeding volatility. Various GARCH family models are found very popular and significant 
to analysis the volatility of stock prices. All these models are extensively applied by the researchers in stock markets and 


still the search of further models is going on. 
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